Dissipation due to tunneling two-level systems 
in gold nanomechanical resonators 
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We present measurements of the dissipation and frequency shift in nanomechanical gold resonators 
at temperatures down to 10 mK. The resonators were fabricated as doubly-clamped beams above a 
GaAs substrate and actuated magnetomotively. Measurements on beams with frequencies 7.95 MHz 
and 3.87 MHz revealed that from 30 mK to 500 mK the dissipation increases with temperature as 
T 0,5 , with saturation occurring at higher temperatures. The relative frequency shift of the resonators 
increases logarithmically with temperature up to at least 400 mK. Similarities with the behavior 
of bulk amorphous solids suggest that the dissipation in our resonators is dominated by two-level 
systems. 

PACS numbers: 85.85.+j, 62.25.-g, 65.40.De 



Nanomechanical resonators with frequencies ranging 
from a few MHz up to 1 GHz have important poten- 
tial applications as ultra-sensitive force sensorsji and are 
also ideal systems in which to study mesoscopic physi- 
cal phenomena involving mechanical degrees of freedom. 
Coupling nanomechanical resonators electrostatically to 
mesoscopic circuit elements^ allows displacement de- 
tection with almost quantum-limited sensitivity to be 
achieved. Furthermore, when cooled to sufficiently low 
temperatures, nanomechanical resonators are expected 
to display quantum coherence^ A detailed understand- 
ing of the processes responsible for dissipation in nanome- 
chanical resonators is crucial if the full potential of such 
devices is to be realized in device applications. Moreover, 
an understanding of how nanomechanical resonators dis- 
sipate energy at low temperatures is an important prereq- 
uisite for developing appropriate models of environmental 
decoherence in such systems^ 

Despite significant experimental^!^ and 
theoretical 6 -^^^^^ efforts a full picture of the 
mechanisms that control the Q-factor of nanomechanical 
resonators has yet to be established, especially at 
low temperatures. Contributions to the mechanical 
dissipation will always come from clamping losses via 
the support s 11 ' 12 and thermoelastic damping,— but 
for beam and cantilever resonators with large aspect 
ratios and at low temperatures, neither is expected to 
be the dominant factor. Rather, defects within the 
crystalline structure of the resonator, or on its surface, 
are expected to play an important role. 1 The low 
temperature dissipation in larger acoustic oscillators, 
especially those made from amorphous solids, has been 
studied extensively and is widely interpreted in terms 
of defects modeled as tunneling two- level systems. 15 
An important prediction of the standard tunneling 
model for amorphous solids is a characteristic power law 
behavior for the dissipation at very low temperatures, 
with a cubic dependence expected for dielectrics and a 
linear power for metals.~ In contrast, only a handful 
of studies of dissipation in nanomechanical resonators 
below 1 K have been conducted*^ However, recent 



low-temperature measurements^ on beam resonators 
with frequencies in the range 5-45 MHz fabricated from 
Si, GaAs and diamond have suggested a T 1 / 3 power law 
dependence for the dissipation. 

Here we present the results of a study of dissipation 
in polycrystalline gold nanomechanical resonators cooled 
down to 10 mK. Metal resonators are being used increas- 
ingly in low temperature experiments since the fabrica- 
tion is simplified when an underlying dielectric layer is 
not required 3 - ! 16 ' 17 and the Q-factor of these resonators 
has been found to be relatively large, although the tem- 
perature dependence of the dissipation has not been stud- 
ied systematically. As metal resonators are readily fab- 
ricated and actuated as monolithic structures, the inter- 
pretation of dissipation measurements is simpler than for 
dielectric ones where metallization is typically required 
to facilitate actuationiSiS A short account of some of 
our results for one of the gold resonators has been given 
elsewhere J£ 

We fabricated nanomechanical resonators consisting of 
doubly-clamped gold beams using a combination of e- 
beam lithography and wet etching, 1 - Polycrystalline gold 
was deposited by thermal evaporation onto a 3 nm thick 
adhesion layer of titanium and the beam released by wet 
etching the underlying GaAs substrate. Our results were 
obtained using two beams: one with length I — 7.5 /im, 
width w = 300 nm, thickness t = 80 nm and a second 
sample with I = 10.5 /im, w — 250 nm, t — 65 nm. A 
SEM micrograph of the longer gold beam resonator is 
shown in Fig. la with a schematic diagram of the mea- 
surement set-up. Using the SEM we estimate that the 
grain size of the polycrystalline gold is 20-40 nm. 

The Q-factor and frequency of the fundamental out of 
plane flexural modes of the beams were measured using 
magnetomotive methods^ The spectral response of the 
beam was obtained using a continuous wave (CW) tech- 
nique by applying a magnetic field (parallel to the plane 
of the wafer) and passing a radio-frequency alternating 
current through the beam. The response to the resulting 
Lorentz force was detected by measuring the induced emf 
(using a heterodyne detection scheme) as the frequency 
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FIG. 1: (a) Detection set-up used to perform the magnetomotive 
measurements. A SEM micrograph of one of the samples is also 
shown, (b) Typical ring-down data for the 3.87 MHz device taken 
at 1 T, the smooth lines are exponential fits to the envelope (the 
data sets are displaced for clarity), (c) Q^ 1 versus B 2 for the 3.87 
MHz device at three different temperatures; the lines are linear fits 
to the data. 



of the current was tuned through the mechanical reso- 
nance. The frequency and Q-factor were then extracted 
by fitting the spectrum to a standard Lorentzian, after a 
linear background was subtracted. Excitation powers of 
-lOOdBm or less were used and care was taken to ensure 
that the resonator remained in the linear response regime 
throughout. 

We also made measurements using a ring-down 
approach 9 in which the resonator was excited by a 200 
/Lts radio-frequency pulses (spaced 1 s apart) with a fre- 
quency close to the mechanical resonance and then the 
response measured in the time domain. Examples of the 
ring-down signal are shown in Fig. lb; the Q-factor was 
extracted by fitting an exponential to the envelope of 
the resulting oscillations. Obtaining accurate estimates 
of the resonator frequency from the ring-down method 
was more problematic because of transient effects which 
mean that the oscillations come closest to the natural 
frequency at later times when the signal itself is rather 
small. 

The resonant frequency of the beams at low tempera- 
tures is affected significantly by differential thermal con- 
traction between the GaAs substrate and the metalj^ii 7 - 
To a good approximation a doubly-clamped beam of 
length I under tension To has a fundamental frequency 
given b y 17 i 20 



fa 



'Et 2 l 2 T Q 
36p + 12ir 2 pA 



(1) 



differential thermal contraction between GaAs and Au 
when the beam is cooled can be estimated usingi 7 - To = 
EA[(^) Au - (^)caAs], where (Al/l) Au is the fractional 
change in length of the gold on cooling. Using Eq. [T]and 
simple estimates of the differential thermal contraction 
obtained from the literatur e) 17 ' 21 we obtain estimates for 
the frequencies of 7.9 MHz and 5.4 MHz for the beams 
with lengths 7.5 and 10.5 pm respectively. The mea- 
sured frequency for the shorter beam was 7.95 MHz, but 
for the longer beam it was 3.87 MHz, somewhat lower 
than our estimate. This deviation is most probably due 
to initial compressive strain which we found (based on 
our experience of fabricating a number of resonators us- 
ing the same process) to be generated in varying degree 
by the fabrication process. 

Use of the magnetomotive transduction technique re- 
sults in a measured Q-factor which is loaded by the exter- 
nal circuitry and the electrical resistance of the sample 
itself. However, the unloaded (intrinsic) Q-factor, Q , 
can be extracted from the magnetic field dependence of 
the measured response, Q££L In our system, the drive 
and detection circuits each have a 50 characteristic 
impedance along with some real resistance due to the 
cables and sample. Using an effective circuit model for 
our experiment in which the resonator is modeled as a 
parallel LCR circuit, we find that (assuming the external 
impedances presented to the sample arc effectively real) 
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where A = wt is the cross-sectional area, E the Young's 
modulus and p the density. The tension arising from 



where Z ext is the external impedance and R m = ''^irfom' 
is the effective resistance of the resonator with £ a mode 
constant of order unity and m the resonator mass. We 
measured the dissipation as a function of B 2 at a range 
of different temperatures. As expected the behavior is 
found to be linear and the gradient allows Qq 1 to be 
obtained by extrapolation to zero field. The results for 
the 3.87 MHz beam are shown in Fig. lc and it is clear 
that the magnetic field can cause a substantial change in 
the dissipation. We found that the gradients showed only 
a very weak temperature dependence (which we ascribe 
to fractional changes in the impedances) . Hence we used 
a simple average of the measured gradients to infer Qq 
at arbitrary temperatures. No corresponding quadratic 
field dependence was found for the resonant frequencies 
of the samples confirming that the impedances seen by 
the sample can be treated as effectively real. 

We carried out a series of measurements of the Q-factor 
and frequency of the 7.95 MHz and 3.87 MHz resonators 
from about 1 K down to 10 mK. The data for the 7.95 
MHz device were obtained from CW measurements at a 
field of 3 T, but for the 3.87 MHz sample CW measure- 
ments were made at 2 T and 3 T, as well as ring-down 
measurements at 1 T. The ring-down measurements at 1 
T provided the most accurate results (as can be seen by 
comparing the scatter in the different data sets in Fig. 3a) 
since less extrapolation was required at this field strength 
to obtain Qq 1 - 
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FIG. 2: 7.95 MHz resonator: (a) Qg as a function of tempera- 
ture, on a log-log scale (main plot) and a linear scale (inset). The 
full line (curve in the inset) is a guide to the eye oc T 0,5 . The 
dashed line is the average of the five points with the highest tem- 
peratures, (b) Fractional frequency shift, Sf/fo, as a function of 
temperature. The line is a logarithmic fit to the data marked with 
full circles (i.e. excluding the two points below 30 mK). Note that 
the data were collected in two runs and a frequency jump of 20.4 
kHz that occurred during thermal cycling between the runs has 
been subtracted from the low temperature data. 



The measured temperature dependence of the intrinsic 
dissipation, Qq , for the resonators is shown in Figs. 2a 
and 3a respectively. The data from the samples follow a 
common pattern and three different regions of behavior 
can be identified. Above about 600 mK the dissipation 
saturates, reaching what appears to be a plateau, whilst 
saturation is also observed for the data below 30 mK. 
In the intermediate region between 30 and 500 mK fits 
to the data show that the behavior of both devices is 
described by a power-law T n , with n — 0.50 ± 0.05 for 
the 7.95 MHz device and n = 0.49±0.02 for the 3.87 MHz 
device. The different sets of measurements of Qq 1 taken 
at different fields and using different techniques overlay 
one another in Fig. 3a. Estimates for the dissipation 
in the plateau regions (shown as dashed horizontal lines 
in Figs. 2a and 3a) are obtained by averaging over the 
five points taken at the highest temperatures for each 
sample. Very similar results are obtained for the two 
devices: = 4.1 X 10~ 5 and Q p \ = 3.8 x 10~ 5 for the 
7.95 MHz and 3.87 MHz resonators respectively. 

The relative shift in the frequency of the resonators 
as a function of temperature is shown in Figs. 2b and 
3b. We have plotted the 2 T CW data for the 3.87 MHz 
resonator as this set provides the best compromise be- 
tween line-width and signal-to-noise ratio. For both sam- 
ples, we find that there is an increase in the frequency 
with temperature over the lower half of the temperature 
range. The behavior here is well described by a logarith- 
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FIG. 3: 3.87 MHz resonator: (a) as a function of temperature, 
on a log-log scale (main plot) and a linear scale (inset). The full 
line (curve in inset) is a guide to the eye <x T 5 . The dashed 
line is the average of the five points with the highest temperatures. 
Only the 1 T ring-down data are shown in the inset for clarity, 
(b) Fractional frequency shift, Sf/fo, as a function of temperature 
(CW data at 2 T). The line is a logarithmic fit to the sub-set of 
data points marked with full circles. 



mic dependence, Sf/fo = C s n)ln(T/To) with gradients 
C s = (7.37±0.95)xl0- 6 andC/ = (8.48±0.30) x 10~ 6 , for 
the 7.95 MHz and 3.87 MHz resonators respectively. The 
frequency shift of the 3.87 MHz device stops increasing at 
about 400 mK and shows a systematic decrease above 600 
mK. However, no decrease is seen for the 7.95 MHz sam- 
ple, though we note that the scatter in the data around 
1 K makes it difficult to draw firm conclusions about 
the underlying behavior at higher temperatures for this 
device. Discrete jumps in the resonator frequency were 
observed when the resonator was warmed up (to about 
4K) between measurement runs. 

The saturation in the dissipation at low temperatures 
could be due to a variety of factors including a crossover 
to a regime where a temperature independent mechanism 
such as clamping losses or Joule heating of the sample 
dominates. However, because some saturation is also ob- 
served in the frequency shift data, a feature seen clearly 
in Fig. 2, we believe that thermal decoupling between 
the sample and the refrigerator is the most likely expla- 
nation. Indeed if we assume the logarithmic dependence 
for the frequency shift in Fig. 2b is obeyed below 30 mK 
and use this as a thermometer for the two points at the 
lowest temperatures then the corresponding dissipation 
points recalibrated in this way lie close to the T 05 line. 

The strong variation seen in the dissipation at such 
low temperatures suggests that tunneling two level sys- 
tems (TLS), which are known to exist in polycrystalline 
metals^ are the dominant source of the dissipation. The 
standard model of low temperature acoustic damping in 
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amorphous solids^ 5 , assumes tunneling TLS with a broad 
distribution of energies and relaxation rates. For res- 
onator frequencies, /, such that hf <C k^,T (as in our 
experiments), dissipation is dominated by the relaxation 
of the TLS (not resonant excitation). A change in the be- 
havior is expected at a temperature T* where the min- 
imum relaxation time for the TLS matches the period 
of the resonator. Above T*, there is a plateau where the 
dissipation ; 15 ' 24 Q~ 1 = nC/2, depends on material prop- 
erties such as the spectral density of TLS, but is indepen- 
dent of temperature and frequency; below T* the dissi- 
pation is expected to have a characteristic power depen- 
dence on temperature that depends both on how the TLS 
relax and the form of the underlying distribution. The 
dissipation we see in the gold resonators fits this pattern 
with both devices showing plateau like features above 
about 600 mK, and a power law dependence ~ T 5 at 
lower temperatures down to about 30 mK. Interestingly, 
a recent calculation^ 4 - predicted a T 5 behavior when the 
relaxation of TLS is due to phonons in a nanomechani- 
cal resonator, though a plateau region is not predicted in 
this case. Furthermore, the dissipation in the plateau re- 
gions is very similar for the two devices although it leads 
to a value C ~ 2.5 x 10^ 5 which is about an order of 
magnitude below the range typical for amorphous solids, 
though similar low values are seen 1 - 5 - in some disordered 
metals as well as in stressed dielectric resonators^ 

The standard tunneling model also predicts a fre- 
quency shift 1 - 5 - that increases logarithmically with tem- 
perature for T < T* (due to the resonant interaction 
between the TLS and the acoustic excitation) with the 
gradient given by C. For T > T* a logarithmic decrease 



in the frequency shift is expected, with a gradient that de- 
pends on the dominant relaxation mechanism. Although 
we see a logarithmic increase in the frequency shift at low 
temperatures, the gradients are only about a third of the 
value of C inferred from the dissipation in the plateau 
regions. At higher temperatures the frequency shifts for 
our two resonators show different behaviors and so firm 
conclusions cannot be drawn. 

In summary, we find that the dissipation follows a weak 
power law dependence Qq 1 oc T ' 5 over the range 30 to 
500 mK and above 600 mK the dissipation saturates. 
The relative frequency shift of the resonators shows a 
logarithmic increase up to about 400 mK. The observa- 
tion of features normally associated with tunneling TLS, 
such as a plateau in the dissipation and a power law 
behavior at lower temperatures together with the loga- 
rithmic change in the frequency shift, suggest that this 
is the right paradigm in this case. Nevertheless, our re- 
sults do not fit with current theories. However, differ- 
ences between our results and theories based on the be- 
havior of bulk amorphous solids are to be expected, not 
just because the rates at which the TLS relax should be 
rather different in nanomechanical resonators^ but also 
because the very small volume of the samples® together 
with the effects of tension^ 5 - may lead to a distribution 
of TLS quite different to that usually assumed for bulk 
solids f^ 5 - 
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